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Abstract 

We introduce a continuum percolation model defined on the points of a d-dimensional ho- 
mogeneous Poisson process. Each Poisson point is connected to all points within its connection 
range, which depends on the distances to the other Poisson points. We show that the new 
model exhibits a phase transition, and obtain results about the critical values in low and high 
dimensions. 
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1 Introduction 

A continuum percolation model consists of a point process and a rule for connecting the points. 
The first mention of these types of models appears to date from 1961 [3j. Typically, the point 
process is a homogeneous Poisson process X on M. d with some given density A, and the points are 
connected in one of two ways. In a random connection model, points are connected to each other 
by undirected edges determined by some probability measure. In a Boolean model, a ball is placed 
around each point in X, and its radius is generated by some probability measure. Two points are 
considered connected if their balls overlap, and clusters are formed in the obvious way. An overview 
of these types of models can be found in [8]. 

In their 1996 paper on continuum percolation [5], Meester and Haggstrom looked at some 
models for which the density of the underlying Poisson process is irrelevant to the percolation 
probability. In these models, altering the density is equivalent to rescaling the model. In this 
paper, we generalise this work. Let us begin by reviewing the original models. 

The hard sphere Boolean models are those Boolean models that assign probability 1 to the 
configurations where none of the balls overlap, although they are allowed to touch. One such 
example is the dynamic lily-pond model, where the configuration of balls is generated by growing 
balls at each point linearly in time. As soon as a ball touches another ball, it stops growing. The 
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existence of this model is non-trivial since there is a.s. no first contact. It has been shown that this 
model exists and that there is no percolation in any dimension [5]. 

The nearest neighbour model is an example of a random connection model. The model is defined 
by the following connection rule: connect each x S X to its k nearest neighbours. Letting U denote 
the event that an infinite component is formed, we define 

k c (d) = mm{k > 1 : F(U) > 0}. 

By ergodicity, we have that k c (d) = min{k > 1 : F(U) = 1}. It has been proven for this model that 
k c (l) = oo and that 2 < k c (d) < oo, for d > 2, see [5]. Furthermore, it is known that k c (d) = 2 for 
large d. It is also believed, but unproven, that k c (2) = 3 and k c (d) = 2 for d > 3. We shall denote 
the nearest neighbour model by NN(d, k), where d is the dimension and k £ N is the parameter. 

In the next section we introduce a generalised nearest neighbour model, of which the NN(d, k) 
model is a special case. The main questions that we shall study in this new model concern the 
percolation probability, and in particular the existence of phase transitions as the dimension of 
the Poisson process (or any other parameter for that matter) is varied. An important tool at 
our disposal is renormalisation. Renormalisation means discretising the space into boxes to form 
a lattice and then looking at the connections between neighbouring boxes. This then permits a 
comparison with site percolation. There are a number of examples of this kind of approach [91 [2]. 
Some of our proofs will be based on the methods used in [5]. 

2 A generalised nearest neighbour model 

Let X be a homogeneous Poisson process on M. d . Consider a point i£l and let di{x) denote the 
Euclidean distance from x to its ith nearest neighbour. The NN(d, k) model can now be defined 
by the connection rule "x is connected to all points within distance d^{x) of x." We want to 
generalise this model by not considering the di(x) themselves, individually, but a function of all 
these distances. Since we want the model to be density- invariant, it is clear that this function 
should be linear in that f({di(x)}i>i) = YliZi fi(d-i(%)), with the {/j}j>i being linear maps from 
M+ to R+, i.e., fi(x) = aix. 

With this motivation, we define a new model by connecting x by an edge to all points y with 
\x — y\ < r(x), where r(x) is given by 

oo 

r(x) = y~]aidi(x), 
i=i 

for some vector a = (a±, oti, ■ ■ ■ ) with on > for all i. We call this model the generalised nearest 
neighbour model, and denote this model by GN(d, a). Since we shall begin our analysis by looking 
at simpler versions of this model, we denote by GNfc(d, a) the model when only the kth component 
of a is non-zero, and has value a. Observe that GN/ C (ci, 1) = NN(d, k). 
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Looking at the GNi(d, a) model, we define a critical value for a as before: 

a c (d) = inf{a > : F{U) > 0}. 

Since the critical value for NN(c£, 1) is known to be oo, we have that a c (d) > 1 for all d. Note 
that an alternative approach would be to fix a and look for the critical value of d. One could also 
consider a corresponding Boolean model, where on each point x of the Poisson process a sphere 
of radius ad\(x) is placed. This model is related to previously considered dependent continuum 
percolation models: for example, for a > 1, the Boolean model dominates all hard sphere models, 
since in a hard sphere model the ball at x must necessarily have radius less than or equal to d\{x). 

For the GN^((i, a) model, we denote the corresponding critical value by ai\d). Hence under 
this notation, a c {d) = ai\d). Sometimes, however, we will vary only k and look for a phase 
transition for GN^(<i, a) with a and d fixed. Note the following property when a < 1: for all points 
x S X, there are at most k — 1 points within distance r(x) = adk(x) < dk{x). This leads to the 
result that the GNfc(d, a) model is dominated by the GNfe_i(d, 1) model for all a < 1. In the 
GN((i, a) model there is no obvious definition of a critical value as we are dealing with a (possibly 
infinite dimensional) vector. However, one number that will be of interest is \a\ = Yli^=i a i- Note 
that percolation is trivial in all dimensions if \a\ = oo. This case is therefore ignored, which restricts 
us to a such that — > as i — > oo. 

In the formation of large clusters, there are two phenomena that are competing. First, to form 
large clusters, points in the cluster should be connected to (many) other points, and to do this, 
points should be close to each other. On the other hand, points should have a large connection 
range, so points should be far away from each other. As we shall see in the next section, this makes 
the model already far from trivial in one dimension. This is in contrast to the NN(d, k) models, for 
which it is very easy to show that there is no percolation in one dimension. 

Although the generalised nearest neighbour model is a random connection model, it is possible 
to define a corresponding Boolean model. In the Boolean version, two points are connected if their 
connection ranges overlap. The proofs of a number of our results can be easily modified to apply 
to this Boolean model, for example Theorem 14.11 below. Typically any numerical bounds on the 
critical value, such as Corollary 14. 2\ will be halved since the connection ranges now combine to 
make connections. 

This Boolean model is of particular interest when we set a\ = 1 and on = for all % > 1. In this 
case, the Boolean version of the generalised nearest neighbour model dominates all homogeneous 
Poisson hard-sphere models. It is an interesting open problem to determine in what dimensions 
percolating hard sphere models exist. The existence of such models has recently been proven for 
W 1 when d > 45 [lj and the existence is also known on some spaces other than W 1 [7]. On M, it 
is known that no such models exist, so that the problem is still open for 2 < d < 44. Showing 
that a c (d) > 1 (for some d) for the Boolean version of the generalised nearest neighbour model 
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would prove that there is no percolation in (i-dimensional hard-sphere models. However, numerical 
simulations for d = 2 suggest that this approach will not be successful. 

The paper is organised as follows. We first study the GN model in one dimension, in Section 
[3l Then in Sections 0] we consider the GNi(<Z, a) model, where only the distance to the nearest 
neighbour is important. In Section [5] we study the GNfc(d, a) model for large d and in Section [6] 
we obtain some results for the general GN(d,a). The paper concludes with the statement of some 
open problems. 

3 The model in dimension 1 

For the nearest neighbour model NN(l,/c) there is no percolation in one dimension for any k. This 
result is a straightforward consequence of the following two features. First, the Poisson process on 
M contains arbitrarily large gaps. Second, if there is no edge between two neighbouring points, then 
these points belong to separate clusters. For the GNfc(l, a) model though, this second property does 
not hold and the proof that there is no percolation in one dimension, given below, is non-trivial. 

We first introduce some notation. Let (x, x+m) be called a gap of length m if Xn(x,x+m) = 0. 
The term m-gap denotes a gap of length greater than m. We say that there is a bridge over a gap 
if two points on different sides of the gap have an edge between them; we say that x bridges a gap 
if the connection range of x, r(x), spans that gap. Furthermore, a gap (x,x + to) is bridged from 
the right if there is a point y > x + to such that y — r{y) < x. The properties of homogeneous 
Poisson processes imply that for any point x £ X, there are infinitely many TO-gaps in both the 
positive (to the right) and negative (to the left) directions for all to < oo. 

Let p(m) denote the probability that an m-gap is not bridged from the right. 

Lemma 3.1 Consider the GNk(l, a) model, and let (3 > (a V 1). Then p((3 2 ) > 0. 

Proof: It will be convenient to use a second distance function, dk(x), where dk(x) is the distance 
from x to its kth nearest neighbour to the right. Furthermore, call x G X a (3-point if df.{x) > (3. 
Observe that if x is not a /3-point, then r(x) < ad~k(x) < a/3 < /3 2 . Hence, when looking for bridges 
from the right over /? 2 -gaps, it is sufficient to only consider (3 points. 

Consider a /5 2 -gap, and the first point to the right of it. Without loss of generality, call this 
point the origin, 0. Let = Xq < X% < X2 < ■ ■ ■ denote all points to the right of and let 
0<Yo<Yi<l2<... denote the /3-points to the right of 0. A sufficient condition for the /3 2 -gap 
to be unbridged from the right is that Yi — ad~k(Yi) > — (3 2 for all i > 0. From the definition of 
/3-points it follows that Y i+ k — Y. L > (3, and hence Yi > J/3- 
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Since the events {adk(Yi) < t} are positively correlated for all i and t, we have 

p(p 2 ) > P{ad k {Yi) <Yi+(3 2 for all i) 

> P{ad k {Yi) < + p 2 for all i) 

oo 

>n p («4(^)< L|:J/? + /? 2 )- (i) 
i=i 

Let r have a gamma ( 1) distribution, and note that d k (x) = r. Since Yi is a /3-point, and 
/3 > (a V 1), we have the following for all i > k(3: 

P(ad k (Yi) < [^\(3 + f3 2 ) > P(d k (Yi) < L^J + 0) > P(d k (Yj < i) 

> P(r < ~ I r > p) = l - p(t > ~ I r > p). (2) 

k k 

We now use that for n — > cxd, 

P(r > n|r > /3) = J ~ c(/3, fc)^"^-", (3) 

where l/c(p,k) = Y,i=ie' f3 (k - Combining (JU), © and ©, we find that p{(3 2 ) > 

n£i(l — a i) with aj ~ c'(a, P, k)i k ~ l e~ l as i — > oo and c'(a,P,k) > 0, and Oj > for all i. Since 
Oj — > exponentially fast, we conclude that p(P 2 ) > 0. □ 

Theorem 3.2 In the GNk(1,o) model, we have a^\l) = oo. 

Proof: We prove the theorem by showing that there exist a.s. infinitely many unbridged /3 2 -gaps 
in both directions. By symmetry, the probability that a /3 2 -gap is bridged from the left is p(P 2 ) as 
well. Now note that for a gap of given length, between two Poisson points x and y say, the event 
that it is unbridged from the right is positively correlated with the event that it is unbridged from 
the left. Indeed, the absence of a bridge from the right makes the distances of points to the right 
of y to their kth nearest neighbour stochastically smaller, and as a consequence, the same holds for 
the distances from the points to the left of the gap to their kth. nearest neighbour. Therefore, the 
probability that a /3 2 -gap is unbridged is at least p(P) 2 , which is strictly positive by Lemma 13.11 

This would complete the proof if /3 2 -gaps were bridged independently. Since this is not the case, 
we describe a scanning procedure to demonstrate that the desired result holds. Assume without 
loss of generality that there is a Poisson point at the origin. Start at the origin and scan to the right 
until you find a /3 2 -gap. This will be unbridged with probability p(P 2 )- If the /3 2 -gap is not bridged, 
then this makes the neighbouring /? 2 -gaps more likely to be unbridged, by the same reasoning as 
above. We then continue scanning to the right of the gap. 

If there is a bridge over the /3 2 -gap, originating at the Poisson point x say, then we jump ahead 
to x + r(x), walk k points to the right, and continue scanning. Since the dependence between the 
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connection ranges of points does not go further than their k nearest neighbours, this removes the 
dependence on the possible bridging of the previous gap(s). 

By the positive correlation (previous gap was unbridged) or independence (previous gap was 
bridged), all gaps encountered by the scanning procedure are unbridged with probability at least 
p{(3) 2 . Hence, we can recover an infinite sequence of /3 2 -gaps such that each is unbridged with 
probability at least p(f3 2 ), and the events that gaps are unbridged are either independent, or posi- 
tively correlated. Hence, there are a.s. infinitely many unbridged /3 2 -gaps to the right. The proof 
is completed by noting that the scanning procedure to the left works exactly the same. □ 

We now show that if the are large enough, then the GN(1, a) model does percolate. In fact, 
we show that for these oti, the model is fully connected. 

Proposition 3.3 In the GNk(l,a) model, let x be a point of the Poisson process on R. If 
= oo, then r(x) = oo a.s. As a consequence, the GN(l,a) model is fully connected. 

Proof: First note that (di(x), d 2 (x) — d\(x),...) = (U\, U 2 , . . .), where Ui,U 2 ,..- are i.i.d. expo- 
nential random variables with parameter 2. Hence, defining do(x) = 0, we may write 

oo oo i oo i 

r ( x ) = S ^oiidi(x) = ^cti^djix) - dj-i(x) = ^ai^Uj 

i=l i=l j=l i=l j=l 

oo oo 

= Ui oh + U 2 Oi + '~ = PiUi + P2U2 + ■■■ , 

i=l i=2 

where (3k = X^i=fc a «- Before continuing, note that YliPi = Yli^ a i- Denote the Laplace transform 
of a random variable X by 4>x- Since the Ui are independent, the Laplace transform of V = ^ PiUi 
satisfies 

OO OO 1 OO Q 

M.) = 11^=111^ = 11 0-5^)- 

i=l i=l i=l 

Since the U\~ are independent and {V = 00} is a tail event, its probability is either or 1, by a 
Zero-One Law. Hence, V = 00 if and only if 4>v( s ) = for all s > 0. But <pv( s ) = if and only if 
fti = 00. Hence, if ^ $ = 00, then V = 00 a.s., and therefore r(x) = 00 a.s. □ 

There is an obvious potential extension of this result to higher dimensions. This result is given 
in Section [6l 

Observe that E(r(x)) = E(E7i) J2Zi a i + E ( C/ 2) ES2 «i + • • • = 3 ESl Therefore, if the 
conditions of Proposition [3731 are not satisfied, the expected connection range is finite. Furthermore, 
it can be shown that the second moment of r(x) can be bounded from above by a polynomial of 
its first moment. Therefore, in this case, the variance is also finite. For related one-dimensional 
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independent percolation models, we have the result that there is no percolation if E(r(x)) < oo [6]. 
Thus, it seems reasonable to conjecture that the converse of Proposition 13.31 is true. 

Conjecture 3.4 IfY^x^i < °°> then all clusters in the GN(l,a) model are finite almost surely. 

Currently, we have only the following weaker Theorem. 

Theorem 3.5 If on = 7* for some < 7 < 1/2, then the GN(l,a) model does not percolate. 

Proof: Since the events that gaps are not bridged are positively correlated, it suffices to show that 
there exists an m such that the probability that an m-gap is not bridged is strictly positive. Since 
the events that a gap is not bridged from the left and from the right are positively correlated, it 
suffices to show that the probability that a gap is not bridged from the right is strictly positive. 

To do that, denote the points of the Poisson process on [0, 00) by Xq < X\ < ■ ■ ■ . Then for all 
i and k we have 

di(X k )<di(X ) + X k -X . 

Hence, we can calculate 

r(X k ) = ^7X(X fe ) < Y, 1%( X ) + X k ~ X ] 

i i 

= r(X ) + -X_ [ Xk - X ] < r(X ) + X k - X . 
1-7 

Hence, for all k, 

X -r(X )<X k -r(X k ). 

So, if Xq does not bridge the gap to its left, then neither do Xi, X2, — To complete the proof, it 
suffices to show that there exists an m such that P(r(Xo) < m) > 0. Letting di denote the distance 
to the ith nearest neighbour on the right, we calculate, 

P(r(X ) <m) = P \J2Ydi(X ) <m^j>P f J^£(X ) < mj 

= P^7?i^'' + P 1 )<™|=P^ <m ) . 
where the Ui are i.i.d. exponentially distributed with parameter 1. Hence, by the Markov inequality, 

i k u k l 7 



(OO k jj \ 1 oo 



7 mil — 7) 2 

We conclude that P{t{Xq) < m) > for m > 7/(1 — 7) 2 , which completes the proof. □ 

Obviously, Theorem 13.51 also holds for all a with < 7* for all i and some < 7 < 1/2. 
Furthermore, readers familiar with long-range percolation may also see superficial similarities with 
the one-dimensional homogeneous case where q(n) < 1 denotes the probability of being connected 
to a vertex at distance n. If ^2 nq(n) < 00, there is no percolation [3]. 
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4 The GNi(d, a) model for d > 2 



Having completed our treatment of the special case d = 1, we continue by looking at d > 2, starting 
with the simplest model: GNi(d, a). The next result shows that there exists a non-trivial critical 
value for all d > 1. 

Theorem 4.1 For d > 2, a c (d) < oo. 

Proof We first consider d = 2, with the density of the Poisson process equal to 1. A 3 x 3 box is 
called a banana box if the lxl box in its centre contains exactly one point, and the rest of box is 
empty, see Figure [TJ The probability that a 3 x 3 box is a banana box is e -1 e -8 . A 3n x 3n box is 
called good if it contains at least one banana box. The probability that a 3n x 3n box is good is at 
least 1 - (1 -e" 9 )" 2 . 

Now choose n so large that the probability that a 3n x 3n box is good is larger than p s ^ te . Con- 
sider two neighbouring good boxes. By construction, in a good box there is a point whose nearest 
neighbour is at distance at least 1. On the other hand, two points in two neighbouring 3n x 3n 
boxes are at most 3ny/5 away from each other. Then for all a > n\/45, any two neighbouring good 
boxes will be connected to each other. Since the probability that a box is good is larger than pf te , 
there is a.s. an infinite cluster of good boxes. Hence, a c (2) < nv45 < oo. For d > 3, the proof is 
similar, and is therefore omitted. □ 



3 3n 
Figure 1: A banana box inside a good box of side length 3n. 

We now use the same idea to find a not so sharp upper bound for a c (2). In the proof of Theorem 
14. II we used unit boxes purely for convenience. To alter the scale of our boxes, we now introduce an 
extra parameter, 5. This 5 allows us to optimise the choice of n, so that we can minimise our upper 
bound for a c (2). A 35 x 35 box is called a <5-banana box if the inner box of side length 5 contains 
one point, and the rest of the box is empty. This happens with probability 5 2 e~ 9s2 . An 35n x 35n 
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box is called 5-good if it contains at least one <5-banana box. So, if we want the probability that a 
35n x 35n box is 5-good to be larger than p c = p^ te , it suffices to have 



we maximise the LHS of dU). Analytically, we readily find that 5 = 1/3. There is also a probabilistic 
way to see this. For a 35 x 35 box to be a 5-banana box, it must contain only one point, and this 
point must be in the central square. Given that there is only one point, the probability that it is 
in the central square is 1/9 for all 5. If we wish to maximise the probability that a 35 x 35 box is a 
^-banana box, we merely need to maximise the probability that there is only one point in the box. 
This leads to 5 = 1/3, so that the smallest n that satisfies H] is given by 



Monte Carlo simulations suggest that the critical value for site percolation on the square lattice Z rf 
is approximately 0.59 [TU]. Since the exact value is unknown, in the following corollary we use the 
best available rigorous upper bound: p c < 0.679492 [12j . 

Corollary 4.2 We have a c (2) < < 41. 

We now move to the asymptotical behaviour of a c (d) for d —* oo. 

Theorem 4.3 We have a c (d) — ► 1 as d — > oo. 

The proof of Theorem l4.3l builds on the work of Meester and Haggstrom [5] and follows a similar 
approach to Section 3 of their paper. They varied k rather than a and looked at the critical value 
of k for the GNfc(ci, 1) model. In their paper, it was proven that k = 2 is the critical value when d 
is large. It should be noted that the result was not originally presented in this form and has been 
expressed here in the terminology of this paper. The technique they use is showing that there is 
percolation in the GN2(d, 1) model for all sufficiently large d. Combining this with an earlier result 
that the GNi(d, 1) model doesn't percolate for any d yields the result. Here, we derive Theorem 
14.31 by showing that for any a > 1 the GNi(d, a) model percolates for d large enough. 

The proof of this result requires a number of steps. Before giving the formal proof, we give 
an outline of the proof and some preliminary results. Our approach will be to construct a point 
process that with non-zero probability gives a subset of the cluster at the origin of GNj (d, a) that 
is infinite. The existence of such a process would implies that the GNi(d, a) model percolates. The 
point process is built up by a sequence of steps. A step is said to be successful if it produces a 
(finite) subset of the cluster at the origin of GNi(d, a) and initiates two subsequent steps. The idea 
is that by taking the dimension very large, we can uniformly bound from below the probability that 



l-(l-5 2 e-^ 2 r 2 >p c . 



(4) 



Since our bound for a c (2) will be n^/45, we want n in ^ to be as small as possible. To do that, 
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a step is successful. By making this lower bound sufficiently large, we then show that an infinite 
sequence of steps occurs with positive probability. 

The out-cluster of a point x in GNi(d, a) is the set of points defined by the following iterative 
procedure, initiated by the set {x}. Given that we have a set {x±,X2, ■ ■ ■ ,x n }, add to this set all 
Poisson points y such that d(xi,y) < ad\(xi) for some i. This is then repeated with the new set 
of points. The procedure either continues forever, in which case the out-cluster is infinite, or stops 
when there are no more points to add. By construction, all points in the out-cluster at x must 
belong to the cluster at x. We shall prove percolation for the GNi(d, a) model by showing that (a 
subset of) the out-cluster at the origin is infinite with positive probability. 

We attempt to construct a subset of the out-cluster at the origin via a point process. This point 
process is built up algorithmically using d-dimensional spatial branching processes. Initially, the 
space (M. d ) is empty and the spatial branching processes place points in this space to create the 
point process. Each spatial branching process (SBP) is run for n generations and is thought of as 
a step in the algorithm. A step of the algorithm is declared successful if two things happen. First, 
the union of the SBP with all previously successful steps must form a subset of the out-cluster at 
the origin of GNi(ci, a). This means that for all points x in the point process, the offspring of that 
point must be within distance ad\{x) of x. Also, the collection of points must obey the law of a 
homogeneous Poisson process. Second, there must be two points in the n-th generation of the SBP 
that satisfy a certain location condition as defined below. These two points are used to initiate two 
more steps of the algorithm. The point process is then defined as the union of all the successful 
steps of the algorithm. 

The SBP lives in R d , but one can consider its projection onto R via the linear map 

L(x 1 , . . . ,x d ) = Vd(x 1 ,x 2 ). (5) 

This permits a comparison with oriented site percolation on the lattice £ = {(i,j) £Z 2 :i> 0, \j\ < 
i, (i + j)/2 £ Z}, with oriented edges from to (i + l,j ± 1). Each site of C corresponds 
to a square S itj = [M(i - \),M(i + \)\ x [M(j - \),M(j + \)\, where M will be chosen later. 
The second condition for a step to be successful (mentioned above) is now that the projected SBP 
starting in box Sij has points in its n-th generation in Sj+ij-i and Sj+ij+i. This allows us to start 
a pair of new SBPs from these points with projected origins in Sj+ij-i and Sj+ij+i respectively. 
In this way, a successful step in the algorithm corresponds to an open site in C Let p c < 1 be the 
critical value of oriented site percolation on C If for all sufficiently high dimensions we can bound 
uniformly (i.e., irrespective of what has happened in previous steps of the algorithm) from below 
the probability of a successful step of the algorithm by some p > p c , , then we have shown that 
the GNi(d, a) model percolates. Before proving this result, we define the spatial branching process 
and give some properties. 
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The spatial branching process. (SBP) Let S r (x) denote the hyper-sphere of radius r centred at x, 
write S r = S r (0), and let \S r \ denote the volume of S r . The spatial branching process with origin 
is defined as follows. We start with Zq = {0}. Given Z n , the offspring of each y G Z n is generated 
by the following procedure. Let 8\ > 0. An independent homogeneous Poisson point process 
on Si+Si{y) 1S generated, with density X(d) such that the expected number of points contained in 
Si(x) is 1. Then a ball is grown around y until either the ball has radius 1 + <5i or C2 points of X™' 
have been encountered, for a certain c% G N. The offspring of y are then the points of X^ that 
are contained in this ball. 

In general, the SBP does not generate a homogeneous Poisson process. Consider two points x 
and y such that d(x, y) < 1 + 5\. It is clear that the overlap Si+fo (x) D Si + Sj (y) is non-empty, and 
that considering the offspring of both x and y together, the density of points is doubled on this 
overlap. However, if we condition on y having no offspring on this overlap, then the joint collection 
of offspring does form a homogeneous Poisson process Si+s^x) U Si+g^y). Thus if the overlap is 
small, then with high probability, namely when y has no offspring on this overlap, we can consider 
the union of the offspring to be a homogeneous Poisson process on S1+81&) U S\ + s 1 (y)- This 
reasoning will be crucial when comparing the SBP to the out-cluster at the origin of the GNi (d, a) 
model. The following standard result shows that in high dimensions, this overlap is negligible. 

Lemma 4.4 If ' x\ and X2 are such that d(x\,X2) > 0.9 and ri, r% G (0.9, 1.1), then 

|5 ri (xi) n s r2 (x 2 )\ _^ 

\S ri (xi)\ 

as d — > oo. 

Observe that for all e > 0, |5i+ £ |/|5i| — > cxd as d — > cxd. The offspring distribution for each 
individual in the SBP is distributed like 7 A C2, where Y is Poisson distributed with parameter 
ISi+jj |/|5i|. Hence, for any < c\ < C2, we can chose 5\ = 6i(d) such that for d sufficiently large, 
E(y A C2) = ci and that 61(d) — > as d — > 00. Note that c\ and C2 do not vary with the dimension. 
Thus with this choice of 5\(d), the offspring of a point x converge weakly (as d — ► 00) to a set of 
points uniformly distributed on the surface of the unit sphere centred at x. If c\ > 1, then the 
SBP is supercritical and the probability of extinction can be made arbitrarily small by taking c\ 
sufficiently large (independent of the dimension). 

The following lemma comes from [11] , and used the map L defined in ([5|). 

Lemma 4.5 Suppose U = (U\, . . . , Ud) is uniform on the surface of Si. Then, as d — > 00, the 
two-dimensional random vector L(U) converges in distribution to the bivariate normal distribution 
N(0, 1) with zero mean and as covariance matrix the identity matrix I. 

The mapping L can be used to map our SBP onto R 2 . Let SBP* denote the limit as d — > cxd of 
the mapped SBP. Since L is continuous, this limit is the same as taking the mapping the limit of the 
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SBP. Thus, by the previous discussion and Lemma 14.51 SBP* is the process where each point has a 
Poisson(ci) distributed number of offspring and these offspring are distributed bivariate normally 
with zero mean and the identity covariance matrix. The same process started at x rather than the 
origin is denoted by SBP*. 

Lemma 4.6 Given e > and c\ sufficiently large, we can find a positive integer Nq and a positive 
number M, such that for all x £ Sij = [M (i - |), M(i + |)] x [M (j - \), M(j + |)], the probability 
that the N$th generation of SBP* X contains at least one point in S^ij-i and at least one point in 
Si + \j + i exceeds 1 — e. 

This result is a small perturbation of a result from [5] and is stated without proof. The process 
considered in [5] is a branching random walk, so it is (slightly different) than the process considered 
here. However, the two processes are very similar and the key difference between the branching ran- 
dom walk and our SBP* is that the branching random walk can never go extinct. The role of c\ in 
Lemma [4.6l is to make the probability of survival for A^o generations sufficiently large, so that the re- 
sult of the lemma applies to SBP* also. We are now ready to prove the main theorem of this section. 

Proof of Theorem 14. 31 Fix a > 1 and consider defining the SBP for a fixed dimension by 
choosing 5\ > as follows. First set 82 > such that 1 — 82 > oT , and then choose <5i such 
that 1 + 5i < (1 — 82)0. Consider running the SBP for n generations and that all points are at 
least distance 1 + 8\ from all others points that aren't its parent. Then, for all x £ SBP, the only 
points in Si + s 1 (x) are the offspring of x coming from the homogeneous Poisson process on S'i + 5 1 (x) 
used to generate them. Thus, in this case, the SBP creates a homogeneous Poisson process on the 
space |J Si+$i (x) , where the union is over all points x in the first n — 1 generations of the SBP. 
Furthermore, consider that for all such x, the offspring are born at least distance 1 — 82 away, i.e., 
all offspring of x appear on the annulus Si+s^x) — Sis 2 ( x )- This implies that d±(x) > 1 — 82 for 
all points x in SBP. Therefore, r{x) > (1 — 82)0. > 1 + 81, for all x in SBP. 

We generate an object that is dominated by the cluster of the site percolation process by running 
the following algorithm. The algorithm consists of the steps (0, 0), (1, —1), (1, 1), (2, —2), . . . , with 
step (i, j) only being carried out if at least one of the steps (i— 1, j — 1) and (£— 1, j + 1) is successful. 
Step consists of a SBP started from a point in S{j and is called successful unless one of the 
following errors occur. 

(a) The spatial branching process fails to reach the two neighbouring boxes in the projected 
space. 

(b) An individual is born within distance 1 — 82 of its parent. 

(c) An individual in the projected SBP is further than Rq from the origin (of the branching 
process). The choice of Rq is given below. 
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(d) An individual is born within distance 1 + <5i of an already generated individual that is 
not its parent. 

A step is stopped if any of the above errors occur. To avoid ambiguity, step is started from 
the point in generation Nq of step (i — l,j — 1) (if successful) or step (i — 1, j + 1) (if successful) 
whose projection is closest to (Mi, Mj). We let Tij denote the cr-algebra generated by the indicator 
functions of the successes of steps (0, 0), (1, —1), (1, 1), . . . , (i, j — 2) of the algorithm. 

This algorithm generates a point process consisting of all the points contained in all the steps. 
This point process is then thinned by removing any points that are born in a previously explored 
part of space. Note that these points are a subset of those that cause type (d) errors and can 
only occur in unsuccessful steps of the algorithm. This thinning procedure ensures that the point 
process is a homogeneous Poisson process on a random subset of M. d . Furthermore, we claim that if 
the probability of a successful step is sufficiently large, this algorithm will with positive probability 
generate a subset of the GNi(<i, a) model that contains an infinite cluster. 

It should be clear from the above discussion that, taken on its own, a successful step generates 
a homogeneous Poisson process on a random subset of M. d . Furthermore, by the law of GNi(d, a), 
the origin of the step is contained in the same cluster as two points that can be used to start two 
subsequent steps. This still holds after the point process is thinned (since the thinning can only 
affect points in unsuccessful steps). Furthermore, subsequent steps can only interfere with an earlier 
successful step with precisely those points that are removed with the thinning. Thus, Theorem 14.31 
follows from the following claim. For all sufficiently large d, and for all £ C, 

P(step (i, j) is successful | ^ij) = p > p c . 

The approach is to show that we can simultaneously make the probabilities of each type of 
error arbitrarily small by suitable parameter choices. The reason for introducing error (c) is to help 
bounding error (d) independently of the history of the algorithm. Let 7 be such that 1 — 47 > p c . 
We proceed by bounding the probability of each type of error from above by 7. 

Lemma 14.61 shows that we can make the probability of error (a) arbitrarily small, for some ci, 
iVo and for all d sufficiently large. We choose c\ and iVo such that P(type (a) error) < 7 for all 
suitably high dimensions. We also fix 1 — 82 > a -1 and c% > c\. Recall from the prior discussion 
about the spatial branching process that for c\ and C2 fixed there exists a 81(d) that defines the 
required branching process. Therefore, we have also fixed 81(d). Furthermore, since 81(d) —* as 
d — > 00, 1 + 8\(d) < (1 — 82)01 in high dimensions, as required. 

Recall that for all e > 0, |Si_ e |/|Si| — ► as d — > 00. Since the maximum number of points in 
a step of the algorithm is now bounded by and 82 is fixed, we make P(type (b) error) smaller 
than 7 by taking the dimension high enough. Next we choose Rq such that the probability that all 
individuals of SBP* are within distance Rq of the origin is at least 1—7- This implies that in high 
enough dimensions, the probability of a type (c) error is also less than 7. 
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It now only remains to bound the probability of a type (d) error. This is the only type of 
error that depends on the history of the process and thus bounding this error is more involved. 
We begin by considering a type (d) error in step caused by an individual from step 

such that M\l (i — i) 2 + (j — j) 2 < To, where To is a constant that will be chosen later on. This 
is straightforward for any To since the number of steps to consider is bounded, meaning that 
the total number of Poisson points is also bounded. Thus, by Lemma 14.41 taking the dimension 
high enough makes the probability of this type of error less than j/2. 

To finish the proof, we consider a type (d) error due to a previous step satisfying 

YM^{i-i) 2 + (j-j) 2 \ =Q>T . 

We define the volume of the step (i, j) to be | U Si + s L (x)\, where the union is over all points x 
contained in the step. Note that the volume of step is uniformly bounded and that the step 
can only fail because of a type (d) error with step if one of its points falls within the 
volume of (i, j). Note that if an error of type (c) occurs, we stop running the step of the algorithm. 
So any point that is born a projected distance greater than Rq from the origin of the step has no 
offspring and is not scanned around. Thus, any point in step that we scan around must be 

at least distance Q — 2Rq from any point in step (i, j). Thus from Lemma 14.51 and the exponential 
decay of the normal distribution, for d large enough, the fraction of the projected volume that falls 
into the circle with radius Rq centred at (Mi,Mj) is less than |S'i+ ( 5 1 |/Q 3 for all large Q. 

Since the number of points such that [M^ (i — i) 2 + (j — j) 2 \ = Q is bounded by a 

constant times Q, and the series Ylq>T 1~~ 2 converges, we can make the total of such volume from 
all points small by choosing To sufficiently large. In particular, we can make this volume so small 
that the probability of each individual in step coming within distance 1 + 6i of another point 
already generated (other than its parent) is less than 7/(2c 2 °). Since there are at most points 
in step we get the desired bound. 

□ 



5 Results for GN k (d, a) for d > 2 

In low dimensions, the difference between GNi and GN^ with k > 1 can be quite pronounced. 
However, in high dimensions this is no longer the case, since there the difference between dk(x) and 
di(x) disappears. This leads to some interesting behaviour. 

First, we discuss an immediate consequence of Theorem 14.31 Recall from Section [2] that for all 
a < 1, the GN2(ci, a) is dominated by the GNi(d, 1) model and that there is no percolation in the 

(2) (2) 

GNi(d,l) model for all d. Hence, a y c '{d) > 1 for all d. Since a c < a c , we have the following 
corollary to Theorem 14.31 

(2) 

Corollary 5.1 We have or c (d) — ► 1 as d — > do. 
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In fact, Theorem 14.31 is expected to hold for the GNfc(d, a) model as well, for all k: 

Conjecture 5.2 For all k > I, a c k \d) — ► 1 as d — > oo. 

The next result shows that for d > 2, any a is sufficient for percolation so long as k is suitably 
large. 

Theorem 5.3 For a > and d > 2 fixed, there exists a k such that the model with = a 
percolates on M. d . 

Before giving the proof, we explain the idea of the proof. We divide M 2 into unit squares to 
create a version of 1? . A square is called good if a certain property, depending only on the Poisson 
process X within the square, holds. In this way we obtain independent site percolation on Z 2 , which 
percolates if the probability that a square is good is greater than p c . It will then be demonstrated 
that if we have percolation of these good squares, then the underlying continuum percolation model 
percolates. Finally, we will demonstrate that we can do this in such a way that the probability of 
a good square is greater than p c . 

Proof: Let a > 0. Let p c (d) denote the critical value of site percolation on 7L d and recall 
that p c (d) < 1 for all d > 2 [JJ. The following argument considers the case when d = 2 and relies 
on p c = p c (2) < 1. The generalisation to higher dimensions is straightforward and is therefore 
omitted. 

We divide each unit square equally into n 2 smaller squares, where n is some odd integer. Let 
X n denote the number of points in a typical subsquare. We call a square good if all its subsquares 
satisfy the following property: 

m 

l<X n <—, 
n z 

where m will be chosen later. Now consider two neighbouring good squares. Without loss of 
generality, let these be [0, l] 2 and [1, 2] x [0, 1]. Define the subsquares 



D; 



n—l ■ n— 1 • | I ri~ 1 n— 1 

2 _|_ _^ 2 _|_ 1 ~r L \ „ f 2 2 



in n n n 



n n 



0<i<n, 



see also Figure [2j By assumption, all the Bi are non-empty, and any point in F>i has not more than 
m neighbours within distance (n — l)/(2n). Hence, for all k > m + 1, all points in a subsquare Bi 
have a connection range of at least a(n — l)/(2n). 

Furthermore, points in neighbouring subsquares are at most distance y5/n apart. By choosing 
n > l+2\/5/a, we ensure that all points in neighbouring subsquares of a good square are connected. 
As a consequence, all points in neighbouring good squares are connected. We can now choose A, 
the density of the Poisson process, to be so large that 

P (x„ = o)<!^. 
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Furthermore, we choose m such that 

n z 2n z 

The above relations imply that the probability that a square is good is at least p c . Since squares 
are good independently of the state of the other squares, and points in good squares are connected, 
the good squares dominate independent site percolation on Z 2 with parameter at least p c . Hence, 
the model percolates. □ 
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Figure 2: Two good squares and their subsquares B\, B2, ■ ■ ■ , -B5 for n = 5. 

Note that Theorem 13.21 Conjecture 15.21 and Theorem 15.31 would together imply that the critical 
value is not monotone in the dimension, d. By Theorem 15 .31 there exists a k such that ac (2) < 0.5. 
Theorem 13 . 2 1 shows that a^\l) = 00 for all k. However, Conjecture [521 would give that a { c\d) -» 1 
as d — ► 00 for all k. In particular, for every k there exists a d > 2 such that ai\d) > 0.5. Thus, 

a{ k \l)>ai k \2)<4 k \d). 

6 Results for GN(d,a) for d > 2 

We now look at the general model, by removing the restriction that a can have only one non-zero 
component. 

Theorem 6.1 Let x be a point of the Poisson process on M. d . Then r(x) = oo a.s. if and only if 

Proof: First choose < e < 1. Let (Ti) be a sequence of disjoint d-dimensional annuli centred at x 
such that for all k the volume of is 1 — e and Uj-LjTj is a ball. Since the volume of a ball of radius 
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r in R d is c(d)r d with c(d) = -K d/2 / {T(d/2 + 1)), the outer radius j k ofT k satisfies c(d)j% = k(l-e), 

Now define Y k to be the number of points in T^, for every k, and let JT& = 1 — Yj.. Then the 

random walk 5 n = Xi H h X n has drift EX k = 1 — .EY* = 1— volume(rfc) = 1 — (1 — e) = e > 0. 

Since such a random walk is transient and converges to oo, there exists an a.s. finite random index 
N such that S n > for all n> N. 

If S n > 0, then Y\ + ■ ■ ■ + Y n < n, and the nth nearest neighbour of x is further away than the 
outer radius j n of the nth annulus, i.e., d n (x) > 7 n = ( ^Tgr" ) 1 ^- Hence, we may write 

iV— 1 j oo 

r ( X ) = X] a *^( X ) - a *^( X ) + ( — Tjf ) ^2 a i il/d - 

i i=l °^ ' i=N 

So, if otii 1 ^ = oo, then r(x) = oo a.s. Analogously, by considering annuli with volume 1 + e, 
we can find the upper bound 

M— 1 -. -, oo 

r{x) < ^di{x) + (-±£) ^ o^ 1 ^, 

i=l ^ ' t=Af 

where M is a.s. finite. Hence, r(x) < oo a.s. if < 00 • This completes the proof. □ 

Since the GN^ model is non-trivial when d > 2, the behaviour of the model (i.e., percolation or 
otherwise) depends on both the tail of a and the individual a^s themselves. For this reason, it is 
clear that there can be no analogy to Conjecture 13.41 in the one dimensional case. 



7 Concluding Remarks 

This paper introduces a new continuum percolation model and proves a number of results about 
this model. In particular, the non-triviality of the model and some features of the critical value 
have been discussed. Obviously, there are still many open problems relating to the model, some of 
which are listed below. 

• The GN(l,a) model does not percolate if the series converges? 

• If l«| = Yli^Li a i = Y^i=i Q i < 1 f° r some k, then does there exist a do such that for all 
d > do the GN(d, a) model does not percolate? This is a generalisation of Conjecture 15. 21 A 
consequence of this result would be the non-monotonicity of the critical value (when varying 
as a function of the dimension) . 

• It might be possible to give a fuller description of the GN(d, a) when \a\ = 1. For example, 
when d is suitably large, the GN2(d, 1) model percolates, but the GNi(d, 1) model does not. 
How does the GN(d, a) model behave when a\ + «2 = 1 and ai = for all % > 3. 
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• This paper has briefly touched upon the related hard sphere models. The author's of this 
paper would like to again highlight the interesting open problem as to whether percolating 
hard-sphere models exist in low dimensions. 

Finally, it would also be interesting to look at other properties of the random graphs created 
by this model than just percolation. 
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